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^^ Abstract Free fall has signed the greatest markings in the history of physics through 

"^ the leaning Pisa tower, the Cambridge apple tree and the Einstein lift. The perspec- 

^ tives offered by the capture of stars by supermassive black holes are to be cherished, 

^I^ because the study of the motion of falling stars will constitute a giant step forward 

Qh in the understanding of gravitation in the regime of strong field. After an account 

on the perception of free fall in ancient times and on the behaviour of a gravitating 
mass in Newtonian physics, this chapter deals with last century debate on the repul- 
sion for a Schwarzschild black hole and mentions the issue of an infalling particle 
C^ velocity at the horizon. Further, black hole perturbations and numerical methods are 

presented, paving the way to the introduction of the self-force and other back-action 
related methods. The impact of the perturbations on the motion of the falling particle 
is computed via the tail, the back- scattered part of the perturbations, or via a radia- 



>^ 

Q^ tive Green function. In the former approach, the self-force acts upon the background 

I— I geodesic; in the latter, the geodesic is conceived in the total (background plus pertur- 
bations) field. Regularisation techniques (mode-sum and Riemann-Hurwitz z func- 

►^ tion) intervene to cancel divergencies coming from the infinitesimal size of the par- 

^-H tide. An account is given on the state of the art, including the last results obtained in 

i-H this most classical problem, together with a perspective encompassing future space 

^O gravitational wave interferometry and head-on particle physics experiments. As free 

^-^ fall is patently non-adiabatic, it requires the most sophisticated techniques for study- 

VO ing the evolution of the motion. In this scenario, the potential of the self-consistent 

^^ approach, by means of which the background geodesic is continuously corrected by 

^^ the self-force contribution, is examined. 
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1 Introduction 

The two-body problem in general relativity remains one of the most interesting 
problems, being still partially unsolved. Specifically, the free fall, one of the eldest 
and classical problems in physics, has characterised the thinking of the most genial 
developments and it is taken as reference to measure our progress in the knowl- 
edge of gravitation. Free fall contains some of the most fundamental questions on 
relativistic motion. The mathematical simplification, given by the reduction to a 
2-dimensional case, and the non-likelihood of an astrophysical head-on collision 
should not throw a shadow on the merits of this problem. Instead, it may be seen as 
an arena where to explore part of the relevant features that occur to general orbits, 
e.g. the coupling between radial and time coordinates. 

Although it is easily argued that radiation reaction has a modest impact on radial 
fall due to the feebleness of cumulative effects (anyhow, in case of high or even 
relativistic - a fraction of c - initial velocity of the falling particle, it is reasonable 
to suppose a non-modest impact on the waveform and possibly the existence of a 
signature), it would be presumptuous to consider free fall simpler than circular or- 
bits, or even elliptic orbits if in the latter adiabaticity may be evoked. Adiabaticity 
has been variously defined in the literature, but on the common ground of the sec- 
ular effects of radiation reaction occurring on a longer time scale than the orbital 
period. One definition refers to the particle moving anyhow, although radiating, on 
the background geodesic (local small deviations approximation), of obviously no- 
interest herein; another, currently debated for bound orbits, to the secular changes in 
the orbital motion being stemmed solely by the dissipative effects (radiative approx- 
imation); the third to the radiation reaction time scale being much longer than the 
orbital period (secular approximation), which is a rephrasing of the basic assump- 
tion. 

But in radial fall such an orbital period doesn't exist. And as the particle falls 
in, the problem becomes more and more complex. In curved spacetime, at any time 
the emitted radiation may backscatter off the spacetime curvature, and interact back 
with the particle later on. Therefore, the instantaneous conservation of energy is 
not applicable and the momentary self-force acting on the particle depends on the 
particle's entire history. There is an escape route, though, for periodic motion. But 
energy-momentum balance can't be evoked in radial fall, lacking the opportunity of 
any adiabatic averaging. The particle reaction to its radiation has thus to be com- 
puted and implemented immediately to determine the effects on the subsequent mo- 
tion. It is a no-compromise analysis, without shortcuts. Thus, the computation and 
the application of the back-action all along the trajectory and the continuous cor- 
rection of the background geodesic, it is the only semi-analytic way to determine 
motion in non-adiabatic cases. And once this self-consistent approach shall be mas- 
tered for radial infall, where simplification occurs for the two-dimensional nature of 
the problem, it shall be applicable to generic orbits. 

It is worth reminding that the non-adiabatic gravitational waveforms are one of 
the original aims of the self-force community, since they express i) the physics closer 
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to the black hole horizon; ii) the most complex trajectories; iii) the most tantalising 
theoretical questions. 

The head-on collisions of black holes and the associated radiation reaction were 
evoked recently in the context of particle accelerators and thereby showing the rich- 
ness of the applicability of the radial trajectory also beyond the astrophysical realm. 
As gravity is claimed by some authors to be the dominant force in the transplanckian 
region, the use of general relativity is adopted for their analysis. 

This chapter reviews the problem of free fall of a small mass into a large one, 
from the beginning of science, whatever this may mean, to the application of the 
self- force and of a concurring approach, in the last fourteen years. There is no pre- 
tension of exhaustiveness and, furthermore, justifiably or not, from this review some 
topics have been disregarded, namely: any orbit different from radial fall; radiation 
reaction in electromagnetism; but also the head-on of comparable masses and Kerr 
geometry, post-Newtonian (pN) and effective one-body (EOB) methods; quantum 
corrections to motion. 

Herein, the terms of self-force and radiation reaction are used rather loosely, 
though the latter does not include non-radiative modes. Thus, the self-force de- 
scribes any of the effects upon an object's motion which are proportional to its 
own mass. Nevertheless, to the term self-force is often associated a specific method 
and it is preferable to adopt the term back-action whenever such association is not 
meant. 

Geometric units (G = c = 1) and the convention (-,+,+,+) are adopted, unless 
otherwise stated. The full metric is given by ga^{t^r) = gaf^i^) -^hg^^(t^r) where 
g(^^ is the background metric of a black hole of mass M and h(^^ is the perturbation 
caused by a test particle of mass m. 



2 The historical heritage 

The analysis of the problem of motion certainly did not start with a refereed publi- 
cation and it is arduous to identify individual contributions. Therefore, an arbitrary 
and convenient choice has led to select only renowned names. 

Aristoteles in the fourth century B.C. analysed motion qualitatively rather than 
quantitatively, but he was certainly more geared to a physical language than his pre- 
decessors. His views are scattered through his works, though mainly exposed in the 
Corpus Aristotelicum, collection of the works of Aristoteles, that has survived from 
antiquity through Medieval manuscript transmission |11|. He held that there are 
two kinds of motion for inanimate matter, natural and unnatural. Unnatural motion 
is when something is being pushed: in this case the speed of motion is proportional 
to the force of the push. Natural motion is when something is seeking its natural 
place in the universe, such as a stone falling, or fire rising. For the natural motion 
of objects falling to Earth, Aristoteles asserted that the speed of fall is proportional 
to the weight, and inversely proportional to the density of the medium the body is 
falling through. He added, though, that there is some acceleration as the body ap- 
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preaches more closely its own element; the body increases its weight and speeds up. 
The more tenuous a medium is, the faster the motion. If an object is moving in void, 
Aristoteles believed that it would be moving infinitely fast. 

After two centuries, Hipparkhos said, through Simplikios | |184| , that bodies 
falling from high do experience a restraining factor which accounts for the slower 
movement at the start of the fall. 

Gravitation was a domain of concern in the flourishing Islamic world between the 
ninth and the thirteenth century by ibn Shakir, al-Blrunl, al-Haytham, al-Khazini. It 
is doubtful whether gravitation was in their minds in the form of a mutual attraction 
of all existing bodies, but the debate acquired significant depth, although not bene- 
fiting of any experimental input. Conversely, in Islamic countries experiments were 
performed as deemed necessary for the development of science. In this sense, there 
was a large paradigmic shift with respect to Greek philosophers, more oriented to 
abstract speculations. 

Leonardo da VincQ stated that each object doesn't move by its own, and when it 
moves, it moves under an unequal weight (for a higher cause); and when the wish 
of the first engine stops, immediately the second stops | |12Q| . Further, in the context 
of 15^^ century gravitation, da Vinci compared planets to magnets for their mutual 
attraction. 

Perception of the beginning of modern science in the early seventeenth century 
is connected on one hand to a popular legend, according to which Galilei dropped 
balls of various densities from the Tower of Pisa, and found that lighter and heav- 
ier ones fell at the same speed (in fact, he did quantitative experiments with balls 
rolling down an inclined plane, a form of falling that is slow enough to be measured 
without advanced instruments); on the other hand, modern science developed when 
the natural philosophers abandoned the search for a cause of the motion, in favour 
of the search for a law describing such motion. The law of fall, stating that distances 
from rest are as the squares of the elapsed times, appeared already in 1604 [SSJ and 
further developed in two famous essays |[86][87| . 

After another century, another legend is connected to Newton |21H who indeed 
himself told that he was inspired to formulate his theory of gravitation | [155| by 
watching the fall of an apple from a tree as reported by W. Stuckeley and J. Conduit. 
The fatherhood of the inverse square law, though, was claimed by R. Hooke and it 
can be traced back even further in the historyof physics. 

The pre-Galilean physics, see Drake |62|^ had an insight on phenomena that is 
not to be dismissed at once. For instance, the widespread belief that fall is unaffected 
by the mass of the falling object shall be examined throughout the chapter, through 
the concept of Newtonian back-action and through the general relativistic analysis 
of the capture of stars by supermassive black holes. 



^ Leonardo spent his final years at Amboise, nowadays part of the French Region Centre, under 
invitation of Frangois I, King of France and Duke of Orleans. 

^ His book presents the contributions by several less known researchers in the flow of time, being 
a well argued and historical - but rather uncritical - account. An other limitation is the neglect of 
non- Western contributions to the development of physics. 
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3 Uniqueness of acceleration and the Newtonian back-action 

One of the most mysterious and sacred laws in general relativity is the equivalence 
principle (EP). Confronted with "the happiest thought" of Einstein's life, it is a relief 
for those who adventure into its questioning to find out that notable relativists share 
this humble opinioqj This principle is variously defined and here below some most 
popular versions are listed: 

I All bodies equally accelerate under inertial or gravitational forces. 
II All bodies equally accelerate independently from their internal composition. 

In general relativity, the language style gets more sophisticated: 

III At every spacetime point of an arbitrary gravitational field, it is possible to 
choose a locally inertial coordinate system such that the laws of nature take 
the same form as in an unaccelerated coordinate system. The laws of nature 
concerned might be all laws (strong EP), or solely those dealing with inertial 
motion (weak EP) or all laws but those dealing with inertial motion (semi- 
strong EP). 

IV A freely moving particle follows a geodesic of spacetime. 

It is evident that both conceptually and experimentally, the above different state- 
ments are not necessarily equivalenj^ although they can be connected to each other 
(for instance: the EP states that the ratio of gravitational mass to inertial mass is 
identical for all bodies and convenience suggests that this ratio is posed equal to 
unity). In this chaper, only the fourth definition will be dealt witljjand interestingly. 



^ Indeed, it has been stated by Synge |194| "...Perhaps they speak of the principle of equivalence. 
If so, it is my turn to have a blank mind, for I have never been able to understand this principle...". 
^ For a review on experimental status of these fundamental laws, see Will's classical references 
(209 210], or else Lammerzahl's alternative view 1 115], while the relation to energy conservation 
is analysed by Haugan |98 ] 
^ For the first definition, it is worth mentioning the following observation 1 194 "...Does it mean 



that the effects of a gravitational field are indistinguishable from the effects of an observer's ac- 
celeration? If so, it is false. In Einstein's theory, either there is a gravitational field or there is 
none, according as the Riemann tensor does not or does vanish. This is an absolute property; it has 
nothing to do with any observer's world-line. Space-time is either flat or curved..." Patently, the 
converse is also far reaching: if an inertial acceleration was strictly equivalent to one produced by 
a gravitational field, curvature would be then associated to inertial accelerations. Rohrlich 1 172j 
stresses that the gravitational field must be static and homogeneous and thus in absence of tidal 
forces. But no such a gravitational field exists or even may be conceived! Furthermore, the particle 
internal structure has to be neglected. 

The second definition is under scrutiny by numerous experimental tests compelled by modem 
theories as pointed out by Damour |44| and Fayet |78 1. 

First and last two definitions are correct in the limit of a point mass. An interesting discussion 
is offered by Ciufolini and Wheeler 1 37 ] on the non- applicability of the concept of a locally inertial 
frame (indeed a spherical drop of liquid in a gravity field would be deformed by tidal forces after 
some time as a state of the art gradiometer may reach sensitivities such to detect the tidal forces of 
a weak gravitational field in a freely falling cabin). Mathematically, locality, for which the metric 
tensor g^y reduces to the Minkowski metric and the first derivatives of the metric tensor are zero, is 
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it can be reformulated, see Detweiler and Whiting |[58 2Q7| , in terms of geodesic 



motion in the perturbed field. Then, the back-action results into the geodesic mo- 
tion of the particle in the metric g(^p (r) + h^n (r, r) where h^n is the regular part of 
the perturbation caused by a test particle of mass m. Thus, the concept of geodesic 
motion is adapted to include the influence of m through h^n . 

A teasing paradox concerning radiation has been conceived relative to a charge 
located in an Earth orbiting spacecraft. Circularly moving charges do radiate, but 
relative to the freely falling space cabin the charge is at rest and thus not radiating. 
Ehlers 1 170] solves the paradox by proposing that "It is necessary to restrict the class 
of experiments covered by the EP to those that are isolated from bodies of fields 
outside the cabin". The transfer of this paradox to the gravitational case, including 
the case of radial fall, is immediate. 

The EP is receptive of another criticism directed at the relation between the foun- 
dations of relativity and their implementation: it is somehow confined to the intro- 
duction of general relativity, while, for the development of the theory, a student of 
general relativity may be rather unaware of i|j 

A popular but wrong interpretation of the EP states that all bodies fall with the 
same acceleration independently from the value of their mass (sometimes referred 
as the uniqueness of acceleration). This view is portrayed or vaguely referred to 
in some undergraduate textbooks, and anyhow largely present in various websites. 
Concerning the uniqueness of acceleration, non-radiative relativistic modes in a cir- 
cular orbit were analysed by Detweiler and Poisson |57|, who showed how the low 
multipole contributions to the gravitational self-acceleration may produce physical 
effects, within gauge arbitrariness (/ = determines a mass shift, / = 1 a centre of 
mass shift). In (57}, the stage is set by a discussion on the gravitational self- force in 
Newtonian theory for a circular orbit. Herein, exactly the same pedagogical demon- 
stration of theirs is applied to free fall. 

A small particle of mass m is in the gravitational field of a much larger mass M. 
The origin of the coordinate system coincides with the centre of mass. The positions 
of M, m and a field point P are given hy ~^,K and "^, respectively (the absolute 
value of "^ is r and mK -\-M~^ = 0). In case of the sole presence of M, the potential 
and the acceleration at P are given by: 

M M 

Mr) = , ^oM = -V^o(r) = -^l^. (1) 

If m is also present, M is displaced from the origin and the potential is: 



limited by the non- vanishing of the Riemann curvature tensor, as in general certain combinations of 
the second derivatives of g^y cannot be removed. Pragmatically, it may be concluded that violating 
effects on the EP may be negligible in a sufficiently small spacetime region, close to a given event. 
^ Again, this opinion is comforted 1 194] "...the principle of equivalence performed the essential 
office of midwife at the birth of general relativity...! suggest that the midwife be now buried with 
appropriate honours...". 
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Since m<^M, eq. ^ is rewritten in the form of a small variation, that is (P{r) 
0o(^) + 50(r) or else 50{r) = 0(r) — <J>o(^); thus: 



50(r) = 



m 



M 



M 
r 



(3) 



^R 



The potential 50(r) determines a field that exerts a force on m, that is the back- 
action of the particle. The singular term ^s diverges, but isotropically around the 
particle position and thus not contributing to the particle motion. Instead, the re- 
maining regular part acts on the particle. Since m <C M, the regular parts of the 
potential and of the acceleration, being V = 5^("^/r), are: 
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At the particle position, the two components of the acceleration are: 
and finally the total acceleration is given by (in vector and scalar form): 
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The Newtonian back-action of a particle of mass m falling into a much larger mass M 
is expressed as a correction to the classical value. This result is more easily derived, 
if the partition between singular and regular parts and the vectorial notation is left 
aside. The force exerted on m is (the origin of coordinate system is made coincident 
with the centre of mass for simplicity, so that mR = Mp): 



and thus 
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(9) 
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It appears from the preceding computations that the falHng mass is slowed down by 
a factor (1 or 2) proportional to its own mass and dependent upon the measurement 
approach adopted. It may be argued that the m/M term arises because the computa- 
tion is referred to the centre of mass, but to shift the centre of mass to the centre of 
M is equivalent to deny the influence of m. 

But instead, what about the popular belief that heavier objects fall faster? Let 
us consider the mass m at height h from the soil and the Earth radius R^\ since 
p + 7? = /z + 7?0, it is found that: 

The mass is now falling faster thanks to a different observer system and the popular 
belief appears being confirmed. On the other hand, in a coordinate system whose 
origin is coincident and comoving with the centre of mass of the larger body M, any 
back-action effect disappears. For J, the distance between the two bodies, it is well 
known that: 

ma = —;^. (13) 

Nevertheless, the translational speed of the moving centre of mass of M (if the latter 
is fixed, any influence of m is automatically ruled out) is depending upon the value 
of m and the same applies to eq. ([9]): it is not possible to find an universal refer- 
ence frame in which the centre of mass moves equally for all various falling masses. 
Thus, the uniqueness of acceleration is result of an approximation, although often 
portrayed as an exact statemenjj or else consequence of gauge choicq^ The unique- 
ness of acceleration holds as long as the values of the masses of the falling bodies 
are negligible. Correctly stated, the principle hardly sounds like a principle: all bod- 
ies fall with the same acceleration independently from their mass, if ... we neglect 
their mass. Although the preceding is elementary, misconceptions tend to persist in 
colleges and higher education. 



^ The difference between fall in vacuum and in the air has been the subject of a polemics between 
the former French Minister of Higher Education and Research Claude Allegre and the Physics No- 
bel Prize Georges Charpak, solicited by the satirical weekly 'Le Canard Enchaine' 1 1 19 1. The Min- 
ister affirmed on French television in 1999 "Pick a student, ask him a simple question in physics: 
take a petanque and a tennis ball, release them; which one arrives first? The student would tell you: 
"the petanque". Hey no, they arrive together; and it is a fundamental problem, for which 2000 years 
were necessary to understand it. These are the basis that everyone should know." The humourists 
wisecracked that the presence of air would indeed prove the student being right and tested their 
claim by means of filled and empty plastic water bottles being released from the second floor of 
their editorial offices ...and asked the Nobel winner to compute the difference due to the air, whose 
influence was denied by the Minister. But in this polemics, no one drew the attention to the New- 
tonian back-action, also during the polemics revamped in 2003 by Allegre 1 1 1 who compared this 
time a heavy object and a paper ball. Such forgetfulness or misconception is best represented by 
the Apollo 15 display of the simultaneous fall of a feather and a hammer |4 1. 
^ During the Bloomington 2009 Capra meeting, this state of affairs was presented as 'the confusion 
gauge'. 
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It is concluded that the Newtonian back-action manifests itself with different nu- 
merical factors possibly carrying opposite sign (from the Pisa tower - 100 pisan 
arms tall - for an observer situated at its feet, Newtonian back-action shows roughly 
as proportional to 1.7 x 10~^^ m/s^ for each falling kilogram). This feature corre- 
sponds to the gauge freedom in general relativity. 

For the latter, when considering perturbations, the energy radiated through gravi- 
tational waves is proportional to n? jM and thus the energy leaking from the nominal 
motion. Therefore, the concept of uniqueness of acceleration is further affected, as 
it will be shown further. 

Finally, it is quoted |[58| that with only local measurements, the observer has no 
means of distinguishing the perturbations from the background metric. In the next 
section, it is shown that the concept of locality or non-locality of measurements 
associated to free fall, even without taking into account radiation reaction, is far 
from being evident and has fueled a controversy for more than 90 years. 



4 The controversy on the repulsion and on the particle velocity 
at the horizon 



143 



162 



114 



The concept of light being trapped in a star was presented in 1783 by Michell 
in front of the Royal Society audience and later by Laplace | 113[|114| . Preti 

describes the close resemblance between the algebraic formulation of Laplace 

and the concept of a black hole, term coined in 1967 by Wheeler [2061. In the last 
century, the Earth, once the attracting mass of reference, was silently replaced by 
the black hole. But, as many centuries before and after Newtonian gravity were 
necessary to formulate motion on the Earth, it should not be a surprise that it is 
taking more than a century to resolve the same Newtonian questions in the more 
complex Einsteinian general relativity on a black hole. 

The existence and the detectability of gravitational waves, the validity of the 
quadrupole formula are among the notorious debates that have characterised gen- 
eral relativity, as described by Kennefick |110| . But closer to the topic of this chap- 
ter is the, surprisingly since almost endless, controversy on the radial motion in the 
unperturbe d Sch warzschild or properly Schwarzschild-Droste (henceforth SD) met- 
ric J64 65 179|P intertwined with the early debates on the apparent singularity at 



^ Rothman (l73| gives a brief historical account on Droste's independent derivation of the same 
metric published by Schwarzschild, in the same year 1916. Eisenstaedt |75| mentions previous 
attempts by DrosteL63J on the basis of the preliminary versions of general relativity by Einstein 
and Grossmann |72|, later followed by the Einstein's works (general relativity was completed in 
1915 and first systematically presented in 1916 |71 1) and Hilbert's 1 100|. Antoci |2| and Lieb- 
scher 1 3 1 emphasise Hilbert's 1 101 1 and Weyl's |205 1 later derivations of solutions for spherically 
symmetric non-rotating bodies. Incidentally, Ferraris, Francaviglia and Reina point to the contribu- 
tions of Einstein and Grossmann |73 1, Lorentz |125J and obviously Hilbert 1 100] to the variational 
formulation. 
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2M and on the belief of the impenetrability of this singularity due to the infinite 
value of pressiorp^at 9/4 M. 

Most references for analysis of orbital motion, e.g. the first comprehensive anal- 
ysis by Hagihara |97| or the later and popular book by Chandrasekhar |33|, don't 
address this debate, that has invested names of the first rank in the specialised early 
literature. 

An historically oriented essay by Eisenstaedt |76| critically scrutinises the rela- 
tion that relativists have with free fal|^ This section does not have any pretension 
of topical (e.g. photons in free fall are not dealt with) or bibliographical complete- 
ness. The questions posed in this debate concern the radial fall of a particle into a 
SD black hole and may be summarised as: 

• Is there an effect of repulsion such that masses are bounced back from the black 
hole? Or more mildly, does the particle speed, although always inward, reaches 
a maximal value and then slows down? And if so, at which speed or at which 
coordinate radius? 

• Does the particle reaches the speed of light at the horizon? 

The discussion is largely a reflection of coordinate arbitrariness (and unaware- 
ness of its consequences), but the debaters showed sometimes a passionate affection 
to a coordinate frame they considered more suitable for a 'real physical' measure- 
ment than other gauges. Further, ill-defined initial conditions at infinity, inaccurate 
wording (approaching rather than equalling the speed of fight), sometimes tortuous 
reasonings despite the great mathematical simplicity, scarce propension to biblio- 
graphic research with consequent claim of historical findings 1 124], they all con- 
tributed to the duration of this debate. The approach of this section is to cut through 
any tortuous reasoning |159| and show the essence of the debate by means of a clean 
and simple presentation, thereby paying the price of oversimplification. 

Four types of measurements can be envisaged: local measurement of time dT, 
non-local measurement of time dt, local measurement of length dR, non-local mea- 
surement of length dr. Locality is somewhat a loose definition, but it hints at those 
measurements by rules and clocks affected by gravity (of the SD black hole) and 



^^ Earman and Eisenstaedt (68| describe the lack of interest of Einstein for singularities in general 
relativity. The debate at the College de France during Einstein's visit in Paris in 1922 included a 
witty exchange on pression (the Hadamard 'disaster'), see Biezunsk i (2 5|. 

^ ^ The translation of the title and of the introduction to section 5 of |76Jserves best this paragraph 
"The impasse (or have the relativists fear of the free fall?) [..] the problem of the free fall of 
bodies in the frame of [..] the Schwarzschild solution. More than any other, this question gathers 
the optimal conditions of interest, on the technical and epistemological levels, without inducing 
nevertheless a focused concern by the experts. Though, is it necessary to emphasise that it is a 
first class problem to which classical mechanics has always showed great concern ... from Galileo; 
which more is the reference model expressing technically the paradigm of the lift in free fall dear 
to Einstein? The matter is such that the case is the most elementary, most natural, an extremely 
simple problem ...apparently, but which raises extremely delicate questions to which only the less 
conscious relativists believe to reply with answers [..]. Exactly the type of naive question that best 
experts prefer to leave in the shadow, in absence of an answer that has to be patently clear to be 
an answer. Without doubts, it is also the reason for which this question induces a very moderate 
interest among the relativists ..." 
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noted by capital letters T, R, while non-locality hints at measurements by rules and 
clocks not affected by gravity (of the SD black hole) and noted by small letters 
t^^^ Therefore, for determining (velocities and) accelerations, four possible com- 
binations do exist: 

• Unrenormalised acceleration d^ r/dt^ ; 

• Semi-renormalised acceleration d^R/dt^\ 

• Renormalised acceleration d^R/dT^\ 

• Semi-renormalised acceleration d^r/dT^. 

The latter hasn't been proposed in the literature and discussion will be limited to the 
first three types. The former two present repulsion at different conditions, while the 
third one never presents repulsion. 

The first to introduce the idea of gravitational repulsion was Droste |[64]|65j him- 
self. He defines: 



dR 



dr 



2M 
r 



(14) 



which, after integration, Droste called the distance 5 from the horizon. This quantity 
is derived from the SD metric posing dt = 0, delicate operation since the relation 
between proper and coordinate times varies in space as explained by Landau and 
Lif shits 1 112 1; thus it may be accepted only for a static observer (obviously the no- 
tion of static observer raises in itself a series of questions, see e.g. Doughty | [6T| , 
Taylor and Wheeler 1 196].). Through a Lagrangian and the relation of eq. ([14]), for 
radial trajectories Droste derives that the semi-renormalised velocity and accelera- 
tion are given by (A is a constant of motion, equal to unity for a particle falling with 
zero velocity at infinity): 



dR 

dt 




(15) 



d^R 

'di^ 



M 



IM 2{dR/dt) 



2M 
r 



M 



1-2A- 



4AM \ 



2M 



, (16) 



where the constant of motion A is given by: 

-1 



^)- 



{dr/dty 1 



2M 
r 



^^ This definition is not faultless (there is no shield to gravity), but it is the most suitable to describe 
the debate, following Cavalieri and Spinelli I30||31||191| and Thirring 1 1971 . 
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From eq. ([16]), two conditions may be derived for the semi-renormalised acceler- 
ation, for either of which the repulsion (the acceleration is positive) occurs for A = 1 
if r < 4M or else dR/dt > ^jtjl^Jl^TMTr. 

Instead in his thesis |64|, Droste investigated the unrenormalised velocity accel- 
eration and for zero velocity at infinity, they are: 



dr 
dt 




(17) 



d^r _ M 

dt^ r^ 



2M 3 {dr/dtf 
~r 2M~ 



r 



M ( ^ 1M\ ( ^ 6M\ 



r^ V r J \ ^ / 



for which repulsion occurs if, still for a particle falling from infinity with zero initial 
velocity, r < 6M or else dr/dt > 1/^3(1 - 2M/r). 

The impact of the choice of coordinates on generating repulsion was not well 
perceived in the early days of general relativity. Further, many notable authors as 
Hilbert fToTpQ2| , Page p56| , Eddington |[69]|, von Laue |203 1 in the German orig 



nal version of his book, Bauer |24|, de Jans | [5T[|53| although indirectly by referring 
to the German version of 1 ,203 1, arrive independently and largely ignoring the exis- 
tence of Droste 's work, to the same conclusions in semi-renormalised or unrenor- 
malised coordinates. 

The initial condition j^ may astray the particle from being attracted by the grav- 
itating mass. Indeed, Droste | [65| and Page 1 156] refer to particles having velocities 
at infinity equal or larger of \/y/2 for the semi-renormalised coordinates and equal 
or larger of 1/a/3 for the unrenormalised coordinates. These conditions dictate to 
Droste and Page that the particle is constantly slowed down when approaching the 
black hole and therefore impose to gravitation an endless repulsive action. 

In the later French editions of his book, von Laue |203 1 writes the radial geodesic 
in proper time, but it is only in 1936 that Drumaux |66| fully exploits it. Drumaux 
criticises the use of the semi-renormalised velocity and considers eq. ([14]) as defining 
the physical measurement of length dR. Similarly, the relation between coordinate 
and proper times (for dr = 0) provides the physical measurement of time dT\ 



^^ Generally, the setting of the proper initial conditions may be a delicate issue e.g. when associated 
with an initial radiation content expressing the previous history of the motion as it will be later 
discussed; or, in absence of radiation, when an external (sort of third body) mechanism prompting 
the motion to the two body system is to be taken into account. The latter case is represented by the 
thought experiment conceived by Copperstock 1 38 1 aiming to criticise the quadrupole formula. The 
experiment consisted in two fluid balls assumed to be in static equilibrium and held apart by a strut, 
with membranes to contain the fluid, until time ^ = 0. Between ^ = and t = ti, the strut and the 
membranes are dissolved and afterwards the balls fall freely. Due to the static initial conditions, 
there is a clear absence of incident radiation, but the behaviour of the fluid balls in the free fall 
phase depends on how the transition from the equilibrium to the free fall takes place. This initial 
dependence obscured the debate on the quadrupole formula. 
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2M 
dT=dl dt. (19) 

Thereby, Drumaux derives the renormalised velocity and acceleration in proper 
time: 



(fR M r m 

for which no repulsion occurs. This approach is followed by von Rabe |204| , 
Whittaker | |2Q8l , Smivasa Rao (T92J, Zel'dovich and Novikov | |2T4| . Nevertheless, 
McVittie, almost thirty years after Drumaux | |142| , still reaffirms that the particle 
is pushed away by the central body as do Treder | 199|, al so in cooperation with 
Fritze pQQ| , Markley p37l , Arifov |[9lp^, McGruder|l4T|. A discussion on radar 
and Doppler measurements with semi-renormalised measurements was offered by 
Jaffe and Shapiro [106 , 107 J . The controversy seems to be extinguished in the 80s, 
although recent research papers still refer to it, e.g. Kutschera and Zajiczek 1 1 1 1 1. 

For the particle's velocity at the horizon, another, though related, debate has taken 
place in some of the above mentioned references as well as in Landau and Lifshits 



TT2| , Baierlein gg, Janis | [TQ8||TQ9| , Rindler (T69J, Shapiro and TeukolskyjTSOJ, 
Frolov and Novikov f84l, Mitra fTsOl, Crawford and Tereno ^, MuUer |152| the 
last ones being recently published. Whether the velocity is c or less, it is still the 
question posed by these papers. 

The further step forward in the analysis of a freely falling mass into a SD black 
hole has taken place in the period from 1957 to 1997. In these forty yearj^ the 
falling mass finally radiates energy (the radiated gravitational power is proportional 
to the square of the third time derivative of the quadrupole moment which is different 
than zero), but its motion is still unaffected by the radiation emitted. The influence 
of the radiation on the motion of a particle of infinitesimal size wasn't dealt with 
until 1997. 



5 Black hole perturbations 



Perturbations were first dealt with by Regge and Wheeler 1 166 167| , where a SD 



black hole was shown to regain stability after undergoing small vibrations about its 



^^ Free fall has also been studied in other contexts. Synge 1 194 1 undertakes a detailed investigation 
of the problem and shows that, actually, the gravitational field (i.e. the Riemann tensor) plays an 
extremely small role in the phenomenon of free fall and the acceleration of 980 cm/ sec is, in fact, 
due to the curvature of the world line of the tree branch. The apple is accelerated until the stem 
breaks, then the world line of the apple becomes inertial until the ground collides with it. 



14 A. Spallicci 



spherical form, if subjected to a small perturbatioif^ The analysis was carried out 
thanks to the first application to a black hole of the Einstein equation at higher order. 
The SD metric describes the background field g^y on which the perturbations 
hf^iy arise. It is given by: 



ds^ = -fl-^yt^+fl.^\ \r^ + ,2 ^^q2 ^ ^jjj2 Q^^2^ (22) 



Eq. (22) originates from the Einstein field equation in vacuum, consisting in the 
vanishing of the Ricci tensor Rj^y = 0. 

Instead, the Regge-Wheeler equation derives from the vacuum condition, but this 
time posed on the first order variation of the Ricci tensor 57?^v = 0- The generic 
form of the variation of the Ricci tensor was found by Eisenhart f74l and it is given 
by: SRf^iy = — 51/^.^ + ^^uB-y where the tensor SFni, variation of the Christoffel 
symbol (a pseudo-tensor), is: Sr^t = 1/2 g^^ih^y.y^hyy.^^ — h^yy), being the per- 
turbation h^xy = Sgi^y. Replacing the latter in the vanishing variation of the Ricci 
tensor, a system of ten second order differential equations in /z^y was obtained. Ex- 
ploiting spherical symmetry, finally Regge and Wheeler got a vacuum wave equa- 
tion out of the three odd-parity equations giving birth to a field that has grown im- 
mensely from the end of th e 5QJ^ 

Zel'dovich and Novikov |213| first considered the problem of gravitational waves 
emitted by bodies moving in the field of a star, on the basis of the quadrupole for- 
mula, thus at large distances from the horizon, where only a minimal part of the 
radiation is emitted. 

Whilst a less known semi-relativistic work by Ruffini and Wheeler ||177| 178| 



appeared in the transition from the 60s to the 70s, it was the work by Zerilli |216f 
|218| , where the source of perturbations was considered in the form of a radially 
falling particle, that opened the way to study free fall in a fully, although linearised, 
relativistic regime at first order. The Zerilli equation rules even-parity waves in the 
presence of a source, i.e. a freely falling point particle, generating a perturbation 
for which the difference from the SD geometry is small. The energy-momentum 
tensor 7)iy is given by the integral of the world-line of the particle, the integrand 
containing a four-dimensional invariant 5 Dirac distribution for representation of 
the point particle trajectory. The vanishing of the covariant divergence of T^y is 
guaranteed by the world-line being a geodesic in the background SD geometry; 



^^ For a critical assessment of black hole stability, see Dafermos and Rodnianski |43 1. 
^^ A well organised introduction, largely based on works by Friedman |83 1 and Chandrasekhar 
(32) , is presented in the already mentioned book by the latter |33|. Some selected publications 
geared to the finalities of this chapter are to be listed: earlier works by Mathews ^138j , Stachel 
|193| , Vishveshvara |202|; the relation between odd and parity perturbations (34) ; the search for 
a gauge invariant formalism by Martel and Poisson 1 140] complements a recent review on gauge 
invariant non-spherical metric perturbations of the SD black hole spacetimes by Nagar and Rez- 
zolla 1 153 1; a classic reference on multiple expansion of gravitational radiation by Thorne 1 198]; 
the derivation by computer algebra by Cruciani |40 41 1 of the wave equation governing black hole 
perturbations; the numerical hyperboloidal approach by Zenginonglu 1215). 
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in this way, the problem of the linearised theory on flat spacetime (for which the 
particle moves on a geodesic of flat space which determines uniform motion and 
thereby without emission of radiation) is avoided. Finally, the complete description 
of the gravitational waves emitted is given by the symmetric tensor /z^y, function of 

r, 0, (j) midt. 

The formalism can be summarised as follows |216- 218p^ Due to the spherical 
symmetry of the SD field, the linearised field equations for the perturbation /z^y 
are in the form of rotationally invariant operator on /z^y, set equal to the energy- 
momentum tensor also expressed in spherical tensorial harmonics: 



Q[h^v]^T^v[5{zu)]. (23) 

where the 5{zu) Dirac distribution represents the point particle on the unperturbed 
trajectory Zm. 

The rotational invariance is used to separate out the angular variables in the field 
equations. For the spherical symmetry on the 2-dimensional manifold on which t^ r 
are constants under rotation in the , sphere, the ten components of the perturbing 
symmetric tensor transform like three scalars, two vectors and one tensor: 



l^tt •) hfr •) hr 



{hte'^ht(j)),{hre'^hr(j)) 



hee he 
h(i)9 hcb 



In the Regge-Wheeler-Zerilli formalism, the even perturbations (the source term 
for the odd perturbations vanishes for the radial trajectory and given the rotational 
invariance through the azimuthal angle, only the index referring to the polar or lati- 
tude angle survives), going as (— 1)^ are expressed by the following matrix: 



d' 



\ 



2M 
r 

sym 
sym 

sym 



HqY 



HiY 
sym 



sym 



-1 



hoY^e 

H2Y hiYe 

r^[KY^GY^ee] 
sym 



hoY^fj) 

hiY^^ 
r^G (70^ -cot 07^) 



r^ sin^ 6 



K+G 



sin^O 



+coieYe 



\ 



(24) 

where Ho^Hi^H2jho^hi^K^G are functions of (r, r) and the / multipole index isn't 
displayed. After angular dependence separation, the seven functions of (r,r) are 
reduced to four due to a gauge transformation for which G = ho = hi =0, i.e. the 
Regge- Wheeler gauge. 

For a point particle of proper mass m, represented by a Dirac delta distribution, 
the stress-energy tensor is given by: 



ra(5 ^ ^ 



m^^5[r-Zu{tW[n], 



(25) 



^^ Two warnings: the literature on perturbations and numerical methods is rather plagued by edito- 
rial errors (likely herein too...) and different terminologies for the same families of perturbations. 
Even parity waves have been named also polar or electric or magnetic, generating some confusion 
(see the correlation table, tab. II, in |218 1). 
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where Zu{t) is the trajectory in coordinate time and u^ is the 4-velocity. 

Any symmetric covariant tensor can be expanded in spherical harmonics |217| . 
For radial fall it has been shown that only three even source terms do not vanish and 
that two functions of (^, r) become identical. Finally, six equations are left with three 
unknown functions Hq =H2-,K^H\. After considerable manipulation, the following 
wave equation is obtained (for radial fall m = and only the index / survives): 



j2iF/(^r) j2^/(^r) 



-Vi{r)^^l(t,r)=Si{t,r), 



(26) 



where r* 
given by: 



2Mln(r/2M— 1) is the tortoise coordinate; the potential V/(r) is 



V/(r): 



2M\ 2X^{X^iy^6X^Mr^ 



ISAM^r+lSM^ 



r J 



r3(Ar + 3M)2 



being A = l/2(/ — 1)(/ + 2). The source Si{t^ r) includes the derivative of the Dirac 
distribution (denoted 5'), coming from the combination of the /z^y and their deriva- 



tives 



\m 



Si 



r(r-2M)_,^ , ^, 



2{r-2M)K 
"r2(A + l)(Ar + 3M)^ 

r(A + 1) - 3M 3Mu^{r - 2M)^ 



2^0 



r(Ar + 3M) 



S[r-Zu{t)] 



(27) 



for u^ = 1/(1 —2M/zu) being the time component of the 4-velocity and K = 
4mY^(2/ + 1)71. The geodesic in the unperturbed SD metric Zu{t) assumes differ- 
ent forms according to the initial conditionq^ herein, the simplest form is given, 
namely zero velocity at infinity. Then, Zu{t) is the - numerical - inverse function of: 






.2MJ 



mrzu_\V2 
3 \2m) 



2Mln 




(28) 



The coordinate velocity Zu of the particle may be given in terms of its position Zu' 

1/2 



Zu = -[l 



2M 

Zu 



2M 

Zu 



(29) 



^^ There is an editorial error, a numerical coefficient, in the corresponding expressions (2.16) in 
|132| and (2.8) in fl331, which the footnote 1 at page 3 in fl391 doesn't address. 
^^ For a starting point different from infinity or a non-null starting velocity, but not their combina- 
tion, see Lousto and Price |1320134|, Martel and Poisson 11391. 
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The dimension of the wavefunction W is such that the energy is proportional to 
Jq W^ dt . The wavefunction, in the Moncrief form |151| for its gauge invariance, 
is related to the perturbations via: 



m,r) 



A + 1 



, r-2M ( I dK'\ 



K -\- I nn — r — — I 



(30) 



where the Zerilli |I216'| normalisation is used for Wi. For computations, this allows 
the choice of a convenient gauge, like the Regge- Wheeler gauge. The inverse rela- 
tions for the perturbation functions K, H2, Hi are given by Lousto |128[|131| : 

r2(Ar + 3M) ^\ r) '' (A + l)(Ar + 3M)r ' 

(31) 

9M3 +9AMV + 3A2Mr2 + A2(A + I)? ?>M^-XMr^Xr^ 

^^~ r2(Ar + 3M)2 ^ r(Ar + 3M) ''^ 

(r-2M)^,,+ 

K-w^(r-2M)[AV2 + 2AMr-3Mr + 3M2] Ku^{r-2Mf , 

r(A + l)(Ar + 3M)2 ~ (A + l)(Ar + 3M) ' ^ ^ 

Ar^-3MAr-3M\ /c^^z, (Ar + M) ?c^Qz,r(r-2M) , 

^ "^ ''' (r-2M)(Ar + 3M) ' (A + l)(Ar + 3M) (A + l)(Ar + 3M) * 

(33) 

Several works by Davies, Press, Price, Ruffini and Tiom no |l46}J48l|174f[176 l|, but 
also by individual scholars like Chung [35 1, Dymnikova ||67| and the forerunners 
of the Japanese school as Tashiro and Ezawa 1 195 1, Nakamura with Oohara and 
Koijma | ,154J or with Shibata 1 182], appeared in the frequency domain in the 70s and 
fewer later on, analysing especially the amplitude and the spectrum of the radiation 
emitted. Haugan, Petrich, Shapiro and Wasserman | |99|l57||18H modeled the source 
as finite size star of dust. 

For an inf ailing mass from infinity at zero velocity, the energy radiated to infinity 
for all modes |47 1 and the energy absorbed by the black hole |48| for each single 
mode, and for all modea^are given by respectively (beware, in physical units): 



99 99 



££[ = 0.0104^, Ef=0.25'^, Y.^t = \mc\ (34) 

while most of the energy is emitted below the frequency: 



^^ The divergence in summing over all / modes is said to be taken away by considering a finite size 
particle ||48]. 
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f.n = Om^. (35) 

Up to 94% of the energy is radiated between 8M and 2M and 90% of it in the 
quadrupole mode. 

Unfortunately, the analysis in the frequency domain doesn't contribute much to 
the understanding of the particle motion, the limitation having origin in the ab- 
sence of exact solutions. A Fourier anti-transform of an approximate solution, for 
instance valid at high frequencies, does not reveal which effect on the motion has 
the neglect of lower frequencies. Thus, the lack of availability of any time domain 
solution has impeded progress in the comprehension of motion in the perturbative 
two-body problem. Although studies on analytic solutions were attempted through- 
out the years, e.g. Fackerell |77|, Zhdanov [2191, Leaver |116 -118|, Mano, Suzuki 



and Takasugi [136J and Fiziev [81J , they were limited to the homogeneous equation. 



6 Numerical solution 

The breakthrough arrived thanks to a specifically tailored finite differences method. 
It consists of the numerical integration of the inhomogeneous wave equation in time 
domain, proposed by Lousto and Price | |133[|13 41 and based on the mathematical 
formalism of the particle limit approximation developed in 1 132] in the Eddington- 
Finkelstein coordinates |[70|[80|. A parametric analysis of the initial data by Martel 
and Poisson has later appeared |139| . Confirmation of the results, among which the 
waveforms at infinity, is contained in p | . 

The grid cells are separated in two categories, according to whether the cell is 
crossed or not by the particle. The latter category, fig.[l] is then formed by the cells 
for which r ^ Zu{t) and the ^{t^r) evolution is not affected by the source. It is 
then sufficient to integrate each term of the homogeneous wave equation. The wave 
operator allows an exact integration: 



// 

J Jc 



{d^,-df)^dA: 

Cell 



-A[W{t^h,r^)^^{t-h,r'')-W{t,r^^h)-W{t,r^-h)]. (36) 

Instead, the product potential- wavefunction is given by: 



// 

J Jc 



V{r)^dA = 

Cell 

V{r)h^[W{t^h,r'')^W{t-h,r'')^W{t,r''^h)^W{t,r''-h)^(^{h^)]. (37) 

The evolution algorithm defines W at the upper cell corner as computed out of 
the three preceding values: 
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^(r + /z,r*) = -^(r-/z,r*) + [^(^r*+/z) + ^(^,r*-/z)] 



l-yV(r) 



. (38) 



ip(t-\-h,r^) 



ij(t,r* -h) 




i){t,r^-\-h) 



Fig. 1 An empty cell never crossed by the particle worldline. 



r* 



For the cells crossed by the particle, a different integration scheme is imposed, 
fig. [2] The product potential-wavefunction is given by: 



J J Cell 



V{r)'FdA = V{r)x 



[A3^(r + /z,r*) +A2^(r - /z,r*) +A4^(^r* +/z) +Ai^(r,r* -/z) + ^{h^)] , 

(39) 
where Ai,A2,A3^4 are the sub-surfaces of the cell. The integration of the source 
term is given b}p^ 
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, (40) 



where ti corresponds to the time of entry of the particle in the cell and to the time of 
departure from the cell; ZuO is the initial position of the particle; E = ^/l —2M/zuO', 
Sigrii = + 1 if the particle enters the cell on the right, — 1 if on the left; Sigrio = + 1 if 



the particle leaves the cell on the right, — 1 if on the left, fig.p 
algorithm, the value of ^ at the upper cell corner is given b} 
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22l 



^^ There are editorial errors in the corresponding expressions (3.6) in 
^^ There are editorial errors in the corresponding expressions (3.9) in 



133 



133 



and (3.4) in 
and (3.5) in 



1139 



1139 



20 



A. Spallicci 



7(r) 



iit-\-h,r^) 



2pit,r^-\-h) 




i){i,r^-h) 



ij(t-h,r*) 
signl = +l \ 5ign2 = -l 

Fig. 2 There are only three physical cases for the particle crossing the t, r cell. 



-f(f-/!,r*) 



W{t,r*-h,] 
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1 + ^(A2-A3) 

V{r) 



W{t + h,r*) = 

1'it,r*+h) 



1-^(A4+A3) 



(Ai +A3) 



1 ^-l^Aa') // S{t,r)dA. (41) 

J J Cell 



For the value of ^ at ^ = /z, the unavailabiHty of ^ at ^ = — /z, is circumvented by 
using a Taylor expansion of ^(r*, — /z) for the initial conditions ^ = 0: 
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W{r\-h) = ^{r\0)-h^ \t=o . (42) 

The setting of initial conditions constitutes a delicate, technical and largely debated 
issue. Apart from the technical difficulty in the numerical implementation, it suffices 
to state how much it is crucial to match the initial radiation conditions, that represent 
the earlier history of the particle, with its position and velocity. For those ZuO starting 
points which are sufficiently far from the horizon, the errors on the initial conditions 
are fortunately not relevant at later times. 

Another numerical issue is the evaluation of the wavefunction and the perturba- 
tions at the position of the particle, but unfortunately not described in the literature 
and too technical for this book. The wavefunction belongs to the C~^ continuity 
clas^and the values before and after the particle position are computed and com- 
pared to the jump conditions posed on the wavefunction and its derivatives | |17H . 
Further, it is necessary to obtain the third derivatives of the wavefunction to de- 
termine the correction to the geodesic background motion of the particle. Given 
the second order convergence of the above described algorithm, this isn't easily 
achieved without recurring to a fourth order scheme | |130| . 



7 Relativistic radial fall affected by the falling mass 
7.1 The self-force 

It has been addressed in the previous section, that the perturbative two-body prob- 
lem, involving a black hole and a particle with radiation emission, has been tackled 
almost 40 years ago. For computation of radiation reaction, it may be worth recalling 
that before 1997, only pN methods existed in the weak field regime. Indeed, it is only 
slightly more than a decade, that we possess methods | |149|p^4 | for the evaluation 
of the self-forcq^in strong field for point particles, thanks to concurring situations. 
On one hand, theorists progressed in understanding radiation reaction and obtained 
formal prescriptions for its determination and, on the other hand, the appearance 
of requirements from the LISA (Laser Interferometer Space Antenna) project | |123| 
for the detection of captures of stars by supermassive black holes (EMRI, Extreme 
Mass Ratio Inspiral), notoriously affected by radiation reaction. 



^^ The Heaviside or step distribution, like the wavefunction of the Zerilli equation, belongs to the 
C~^ continuity class; the Dirac delta distribution and its derivative belong to the C~^ and C~^ 
continuity class, respectively. 

^^ A point-like mass m moves along a geodesic of the background spacetime if m ^ 0; if not, the 
motion is no longer geodesic. It is sometimes stated that the interaction of the particle with its own 
gravitational field gives rise to the self-force. It should be added, though, that such interaction is 
due to an external factor like a background curved spacetime or a force imposing an acceleration 
on the mass. In other words, a single and unique mass in an otherwise empty universe cannot 
experience any self-force. Conceptually, the self-force is thus a manifestation of non-locality in the 
sense of Mach's inertia 11351. 
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Such factors, theoretical progress and experiment requirements, have pushed the 
researchers to turn their efforts in finding an efficient and clear implementation of 
the theorists prescription^ by tackling the problem in the context of perturbation 
theory, for which the small mass m corrects the geodesic equation of motion on a 
fixed background via a factor ^(m) (for a review, see Poisson p?) and Barack p5|). 

Before the appearance of the self-force equation and of the regularisation meth- 
ods, the main theoretical unsolved problem was represented by the infinities of 
the perturbations at the particle's position. After determination of the perturbations 
through eqs. ( [26l[3Tp3] ), the trajectory of the particle could be corrected simply by 
requiring it to be a geodesic of the total (background plus perturbations) metric (the 
Christoffel connection pji refers to the full metric): 

but the perturbation behaves as: 

h R r^ . . (44) 

'{gyS^UyUs) [xy-zl{T)\ [x^-zf(T)] ' 

thus diverging as the inverse of the distance to the particle and imposing a singular 
behavior to F^n on the trajectory of the particle. Thus, the small perturbations as- 
sumption breaks down near the particle, exactly where the radiation reaction should 
be computed. 

The solution was brought by the self- force equation, formulated in 1997 and 
baptised MiS aTaQuWJ^ from the surname first two initials of its discoverers, who 
determined it using various approaches, all yielding the same formal expression. In 
the MiSaTaQuWa prescription, the self- force is only well defined in the harmonic 
(de Bonder) | [49l[50| gauge (stemmed from the Lorenz gauge 1 126|) and any depar- 
ture from it - its relaxation - undermines the validity of the equation of motion. 

Mino, Sasaki and Tanaka 1 149] used two methods, namely the conservation of 
the total stress-energy tensor and the matched asymptotic expansion. The former 
generalises the analysis of DeWitt and Brehme | [59| and Hobbs | |104| , consisting in 
the calculation of the electromagnetic self-force in curved spacetime previously per- 
formed in flat space by Dirac | ,6QJ. It evaluates the perturbation near the worldline 
using the Hadamard expansion ||96| of the retarded Green function |92j|94|. Then, it 
deduces the equation of motion by imposing the conservation of the rank-two sym- 



^^ It is currently believed that the core of most galaxies host supermassive black holes on which 
stars and other compact objects in the neighbourhood inspiral-down and plunge in. Gravitational 
waves might also be detected when radiated by the Milky Way Sgr*A, the central black hole of 
more than 3 million solar masses |29 82 1. The EMRIs are further characterised by a huge number 
of parameters that, when spanned over a large period, produce a yet unmanageable number of 
templates. Thus, in alternative to matched filtering, other methods based on Covariance or on Time 
& Frequency analysis are investigated. If the signal from a capture is not individually detectable, it 
still may contribute to the statistical background fT7) . 
^^ In 2002 at the Capra Penn State meeting by Eric Poisson. 
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metric total stress-energy tensor, via integration of its divergence over the interior 
of a thin world tube around the particle's worldline. 

The latter, reformulated by Poisson 1 158], in a buffer zone matches asymptoti- 
cally the expansion of the black hole perturbed background by the particle with the 
expansion around the particle distorted by the black hole. 

Also in 1997, the axiomatic approach by Quinn and Wald |1641 was presented. 
To them, the self-force is identified by comparison of the perturbation in curved 
spacetime with the perturbation in flat spacetime. The procedure allows elimination 
of the divergent part and extraction of the finite part of the force. 

On the footsteps of Dirac's definition of radiation reaction, in 2003 Detweiler and 
Whiting [ [58| , see also Poisson | |158| |, offered a novel approach. In flat spacetime, the 
radiative Green function is obtained by subtracting the singular contribution, half- 
advanced plus half-retarded, from the retarded Green function. In curved spacetime, 
and in the gravitational case, the attainment of the radiative Green function passes 
through the inclusion of an additional, purposely built, function. The singular part 
does not exert any force on the particle, upon which only the regular field acts | [54| . 
The latter, solely responsible of the self- force, satisfies the homogeneous wave equa- 
tion and may be considered a radiative field in interaction with the particle. This 
approach emphasises that the motion is a geodesic of the full metric and it implies 
two notable features: the regularity of the radiative field and the avoidance of any 
non-causal behaviouiPI 

Gralla and Wald have attempted a more rigorous way of deriving a gravitational 
self- force |91 1. Their final prescription, namely self-consistency versus the first or- 
der perturbative correction to the geodesic of the background spacetime, shall be 
addressed later in this chapter. On the same track of improving rigour, an alternative 
approach and a new derivation of the self- force have been proposed by Gal'tsov and 
coworkers |88| and by Pound 1 160], respectively. 

The determination of the self-force has allowed not only targeted applications 
geared to more and more complex astrophysical scenarios, but also fundamental 
investigations: on the role of passive, active and inertial mass by Burko |28|; the 
already quoted papers on the Newtonian self-force JSTJ, on the EP 1 5 8, 207], on the 
relation to energy conservation by Quinn and Wald 1 165]; on the relation between 
self- force and radiation reaction examined through gauge dependence and adiabatic- 
ity by Mino |144f[147| ; the differentiation between adiabatic, secular and radiative 



approxi matio ns as well as the relevance of the conservative effects by Pound and 
Poisson |l6l| ; on the relation between h^^^^ and /z^, tail and regular parts of the field 
by Detweiler |[55|. 

The following wishes to be constrained to a physical and sketchy picture of the 
self-force. For this purpose, the original MiSaTaQuWa approach - the force acts on 
the background geodesic - is more intuitive. One pictorial description refers to a 
particle that crosses the curved spacetime and thus generates gravitational waves. 
These waves are partly radiated to infinity (the instantaneous part) and partly scat- 
tered back by the black hole potential (the non-local part), thus forming tails which 



^^ Given the elegance of this classic approach, the self-force expression should be rebaptised as 
MiSaTaQuWa-DeWh. 
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impinge on the particle and give origin to the self- force. Alternatively, the same phe- 
nomenon is described by an interaction particle-black hole generating a field which 
behaves as outgoing radiation in the wave-zone and thereby extracts energy from the 
particle. In the near-zone, the field acts on the particle and determines the self-force 
which impedes the particle to move on the geodesic of the background metric. The 
total force is thus written as: 



Ff"uni^) 



F" 

inst. 



^tail ' 



(45) 



where F^^^ is computed from the contributions that propagate along the past light 
cone and F^^-^ has the contributions from inside the past light cone, product of the 
scattering of perturbations due to the motion of the particle in the curved spacetime 
created by the black holq^ fig.[3](pictorially, in a curved spacetime, the radiation is 
not solely confined to the wave front). The self- force is then computed by taking the 
limit F^^^r = F^^.^ [x^ Zu{t)\. Thus, it is conceived as force acting on the background 



geodesic 1 149 1641, wherein F^ refers to the background metric: 



/A 



Xo,^ X, 



..-•> 



Fig. 3 The radiation is going to infinity (instantaneous) or is scattered back (tail). The latter part 
determines the self-force. The self-force is defined for xio ^ Zu where Zu is the position of the 
particle on the worldline T, while x is the evaluation point. 






d\ 



dx^ 



"" r^.M^uy. 



ipyu 



(46) 



All MiSaTaQuWa-DeWh approaches produce the same equation for the self- 
acceleration, given by: 



ig"^ +u"u^) hsh},- I'^pKs) «'«^ (47) 

where the star indicates the tail (MiSaTaQuWa) or radiative (DeWh) component. 



a%if^-'o«P^.." 



Eq. (47 ) is not gauge invariant and depends upon the de Bonder gauge condition: 



^^ Detweiler and Whiting (58] refer to the contribution inside the light cone via the Hadamard 
expression (96| of the Green function. 
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h^r =^. (48) 

where /z*^ = hZ^ - \gysh^ and /z* = g^^h"^^. Barack and Ori have shown lETl that 
under a coordinate transformation of the form x^ ^ x^ — t,^, under which the per- 
turbation transforms according to: 

Vv -^ Vv + ^m;v + ^v;m , (49) 

the self-force acceleration transforms as: 



■^self ~^ ^self ' 



f^^u^u^)^^R^^^^u^^\'' 



(50) 



where the terms are evaluated at the particle and ^^;lv ^^ ^^^ Riemann tensor of the 
background geometry. Thus, for a given two-body system, the MiSaTaQuWa-DeWh 
acceleration is to be mentioned together with the chosen gaugq^ 

The identification of the tail and instantaneous parts was not accompanied by a 
prescription of the cancellation of divergencies, which indeed arrived three years 
later thanks to the mode-sum method by Barack and coworkers |T3J[l4][T6]|20). The 



mode-sum method relies on solutions to interwoven difficulties, mostly related to 
the divergent nature of the problem, but tentatively presented as separate hereafter. 

Spherical symmetry allows the force to be expanded into spherical harmonics 
and turns out to be once more the key factor for black hole physics, after having 
been the expedient for the determination of the wave equation. The divergent nature 
of the problem is then transformed into a summation problem. For each multipole, 
the full force is finite and the divergence appears only upon infinite summing over 
/,m. 

Furthermore, the tail component can't be calculated directly, but solely as dif- 
ference between the full force and the instantaneous part; thus, the self-force is 
computed as: 



F"lf = }^i:Wlf{x)-F-J^{x) 



X^Zu ^ 



(51) 



Each of the two quantities F?^ (x) and fj-^/^ is discontinuous through the particle 
location and the superscript ± indicates the two (different) values obtained by taking 
the particle limit from outside (x^z'^) and inside (x ^ z~). However, the difference 
in eq. ( [ST] ) does not depend upon the direction from which the limit is taken. 

The full and the instantaneous parts have the same singular behaviour at large 
/ and close to the particle; their difference should be sufficient to ensure a regular 
behaviour at each /. Unfortunately, another obstacle arises from the difficulty of cal- 
culating the instantaneous part mode by mode. Therefore, the divergence is dealt 
with by seeking a function H^^, such that the series J^/ 
gent. 



T^al± _ TTal± 



IS conver- 



^^ The self-force being affected by the gauge choice, the EP allows to find a gauge where the self- 
force disappears. Again, as in Newtonian physics, such gauge will be dependent of the mass m, 
impeding the uniqueness of acceleration. 
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The function H^^ mimics the instantaneous component at large / and close to the 
particle. Once such condition is ensured, eq. ^5^ is rewritten as: 



^self 



ral± 
full 



where 
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'«'±(x) = lim y 
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(52) 



(53) 



The addition and subtraction of the function H^^ guarantees the pristine value of 
the computation. In general, for L = / + 1/2: 



fjal ^^a^_ 



B^ 



-C'/L. 



(54) 



Thus, the mode-sum amounts to 1 14|: i) numerical computation of full modes; 
ii) derivation of the regularisation parameters A^B^C, and D (obtained on a local 
analysis of the Greens function near coincidence, x^Zu,^^ large 1); iii) computation 
of eq. ( 52 ) whose behaviour has to show a 1/L^ fall off if previous steps are correctly 



carried out. 



7.2 The pragmatic approach 



The straightforward pragmatic approach by Lousto, Spallicci and Aoudia | 128||129 
|187|188| is the direct implementation of the geodesic in the full metric (background 
+ perturbations) and it is coupled to the renormalisation by the Riemann-Hurwitz ^ 
function. These two features justify the pragmatic adjective. Though the application 
of the ^ function is somewhat artificial and the pragmatic method is somewhat naive, 
the latter has the merit of: i) a clear identification of the different factors participating 
in the motion; ii) potential applicability to any gauge and to higher orders of the ^ 
function renormalisation. 

Dealing only with time and radial components, two geodesic equations can be 
written and then combined into a single one, after elimination of the geodesic pa- 
rameter. Thus, for radial fall the coordinate acceleration is given by the sole radial 
component: 









re- 



' ht^P ' 



- 2if^Zp - 



-Ftf 



(55) 



where Fq refers to the full metric and Zp is given by eq. ( 56 ). Eq. ( 55 ) refers to 



• the full metric field ^^y (^, r) previously defined; 

• the displacement Az, difference between the perturbed Zp {t) and the unperturbed 
Zu{t) positions, and the coordinate time derivatives: 



■-Zu^Az, 



■-Zu^Az, 



■-Zu^A'z; 



(56) 
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the Taylor expansion of the field and its spatial derivative: 

gliv \zp=8!iv \zu{t) ^^Zgixv.r \r=zu{t)^ 
gjiv^r \zp=gliy,r \zu{t) ^^Zgfxv^rr \r=Zu{t) 



(57) 



The unperturbed trajectory of the particle Zu{t) is given by the inverse of the relations 
r(r), e.g. eq. (28 ). Supposing that the relative strengths of the perturbations and the 



deviations behave as: 



[Mi)]2 
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At 
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Az 
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(58) 



Then, the coordinate acceleration correction is given by a n expansio n up to V^ order 
for all quantities, which corresponds to the expression in |128[|l29p ^ 



Az = ai{g,zu) Az^ a2{g,Zu) Az^ ae(h,zu) 



(59) 



The particle determines in first instance the emission of radiation h(^^ , which af- 
ter backscattering by the black hole potential, interacts with the particle itself result- 
ing into a change in acceleration (the coefficient a^ depending on /z^y and deriva- 
tives). The latter places the particle elsewhere from where it should have been, that 
is Zu{t). The field is thus to be evaluated at this new position resulting into a further 
variation in acceleration (the terms a\Az and a2 At depending on g^y and deriva- 
tives). 



All terms in eqs. (59 62) are of \/M order; the terms a\Az and a2At repre- 
sent the background field evaluated on the perturbed trajectory; a6 represents the 
perturbed field on the background trajectory. The expressions in tab.l are gauge in- 
dependent, while in tab. 2 they are shown in the Regge- Wheeler gauge (Ho = H2 and 
^ = as in head-on geodesies). Finally, the coefficient Oq is the lowest order term 
corresponding to a particle radially falling into the SD black hole and not affected 
by the perturbations. It corresponds to the unrenormalised acceleration and it is to 
be added to the terms of eqs. ( |59|62| ) to compute the total acceleration: 



^ = g^^gttAl - 2^''grrAl + 2^''gtt, 



M 



2M ( IM 
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r \ r 



(60) 

If Az receives its main contribution from the background metric g^y or else cu- 
mulative effects are let to grow, a different expansion may be considerecrM 



^^ Apart from some editorial errors therein, ai,2,6 correspond to the A^B^C coefficients in 128 

11291, which are not to be confused with the A, 5, C coefficients of the mode- sum! 

^^ Supposing that the relative strengths of the perturbations and the deviations behave as: 



[h^ 



1)12 



/z(2) /z(i) Az Az 

< < — ^— , 

8 g Zp Zp 



(61) 
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Table 1 Representation of the terms of eq. ( 59l in gauge independent form. The elements in each 
column are produced by the terms of eq. 1 55}, in the first row. From the algebraic sum of the 
elements between horizontal lines, the terms of eq. ([59} are derived. 
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Table 2 Representation of the terms of eq. (|59| ) in Regge- Wheeler gauge. The elements in each 
column are produced by the terms of eq. ( |55} , in the first row. From the algebraic sum of the 
elements between horizontal lines, the terms of eq. ([59} are derived. 
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then, the coordinate acceleration correction would be given by an expansion up to V^ order in 
perturbations and 2^^ order in deviation 1 188 1: 



Az = ai (g.zu) Az + a2{g,Zu) Az + as (g.zu) Az^ + a4{g) Az^+ 

a5{g,Zu)AzAz^ae{h,Zu)^ccT{h,Zu)Az^a8{h,Zu)At (62) 

In eq. l [62} : i) solely second order terms in perturbations are not considered; ii) the terms 
0C2 (g, Zu) Az, 0C3 (g, Zu) Az^, 0C4 (g) Az^, as (g, Zu) AzAz represent the background field evaluated on 
the perturbed trajectory at second order in deviation; iii) a^s tend to infinity close to the hori- 
zon, conversely to the ai,2 coefficients; iv) a-j {h,Zu) Az, Ofg ih.Zu) ^z represent the perturbed field 
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In radial fall, it has been indicated by two different heuristic arguments | 128||131| | 



that the metric perturbations should be of C^ continuity class at the location of the 
particlq^ One argument 1 128] is based on the integration over r of the Hamiltonian 
constraint, which is the tt component of the Einstein equations (eq.[C7a] in ||218|); 
the other |13H on the structure of selected even perturbations equations. In [^ a 
stringent analysis on the C^ continuity is pursued in terms of the jump conditions 
that the wavefunctions and derivatives have to satisfy to guarantee the continuity of 
the perturbation^^ Anyhow, the connection coefficients and the metric perturbation 
derivatives have a finite jump and they can be computed as the average of their 
values at Zj^ ± e with e ^ 0. 

The supposed C^ continuity class of the metric perturbations allows to deal with 
the divergence with / of the a^ coefficient fl28','1291. The divergence originates 
from the infinite sum over the finite multipole component contributions. One way 
of regularising this sum is to subtract to each mode precisely the / ^ oo contribu- 
tion, since for ever larger / the metric perturbations tend to some finite asymptotic 
behaviour. Thus, the subtraction from each mode of the / -^ oo part leads to a conver- 
gent series. The renormalisation by the Riemann-Hurwitz ^ function was proposed 
first in |128[p^| and then extended to higher orders in |188j|. For L = / + 1/2, it 



can be shown that: 



ae = Y, (^L oci = <±L + a^L^ + a'^^L'^ + a^L'^ + ^{L'^). (64) 



/=o 



Eq. ([64]) is casted to have a similar form to the mode-sum expression. The average of 
a^^ and a^^ vanish at the position of the particle, whereas Y4L0 ^5 ~ °° determines 
the divergence. 

The Riemann ^ function |168| and its generalisation, the Hurwitz ^ function 
| |105| , are defined by: 



on the perturbed trajectory, and the av-g coefficients are larger near the horizon. T hese last two 
coefficients may be regularised in / by the Riemann-Hurwitz ^ functi on as shown in |188j. 
^^ The jump conditions were also dealt with by Sopuerta and Laguna 186 . 
^^ Having suppressed the / index for clarity of notation, after visual inspection of eq. (|26J, con- 
taining a derivative of the Dirac delta distribution, it is evinced that the wavefunction W is of C~^ 
continuity class and thus can be written as: 

W{t,r) = W+(t,r) 0i+W-{t,r) 02, (63) 

where &i = & [r — Zu{t)], and &2 = & lzu{t) — f] are two Heaviside step distributions. Computing 
the first and second, space and time and mixed derivatives, Dirac delta distributions and derivatives 
are obtained of the type 5[r — Zu{t)] and 5'[r — Zu{t)], respectively. It is wished that the disconti- 
nuities of W and its derivatives are such that they are canceled when combined in K, H2 and Hi . 
After replacing W and its derivatives in eqs. |3T|33j, continuity requires that the coefficients of &i 
must be equal to the coefficients of 02, while the coefficients of d and d' must vanish separately. 
After some tedious computing and making use of one of the Dirac delta distribution properties: 
f{r)5'[r-zu{t)] = f[zu{tW[r-Zu{t)]-f'{zu{t))5[r-Zu{t)], at the position of the particle, the 
jump conditions for W and its derivatives are found. Furthermore, the jump conditions allow a new 
method of integration, as shown by Aoudia and Spallicci (61 . 
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CW = I(/r, C(5,a) = £(/ + ar, (65) 

1=1 1=0 

where in our case a = 1/2. Two special values of the Hurwitz functions, namely 
^(0, 1/2) = and ^(2, 1/2) = 1/2 n^, cancel the divergent term and determine that 
the term Y4L0 ^e^~^ g^ts a finite value, respectively. Barack and Lousto |fl8| have 
shown the concordance of the mode- sum and the ^ regularisations for radial fall. 



8 The state of the art 

It is now time to discuss the state of the art of the radial fall affected by its mass and 
the emitted radiation. As shown in the introduction, the adiabatic approximation re- 
quires that a given orbital parameter q changes slowly over time scales comparable 
to the orbital period P (this is somewhat a coarse definition since the small mass 
always 'reacts' immediately): Aq = qP ^ q. For circular and moderately elliptic 
orbits, the above condition, where q is function of the semi-lactus rectum p and ec- 
centricity e, is transformed into a condition on the m/M ratio |42|. In radial fall, 
though, it is far from being evident, and even possible, to identify a condition on 
adiabaticity within which any simplification may occur. The feebleness of cumula- 
tive effects for radiation reaction does not imply their non-existence. On the con- 
trary, this is the case where most care and sophisticated techniques are demanded 
for the computation of the motion affected by the back-action, even if the latter has 
moderate effects. Therefore it is not surprising, thanks to the feebleness and to the 
difficulties, that solely two studies (one based on the pragmatic method, the other 
on the self-force) exist. 



8.1 Trajectory 

The perennial question on the behaviour of the infalling mass reflects itself in the 
determination along which direction the back-action is exerted. 

Lousto (fig. 2a in | |128| |) suggests that the a^ term, denoted therein C (the vari- 
ation of the coordinate acceleration of the particle due solely to perturbations; it 
corresponds to the self-force when referred to coordinate time, see next section), in- 
creases approaching the Zerilli potential at 3.1M and reaches its peak value around 
2AM. The same reference (fig. 2b in |128| ) shows that the coordinate acceleration, 
thus including aiAz and a2 At terms, is slowed dowip] and mostly until around 
3.1M. The two statements are not contradictory if the discrepancy is attributed to 
aiAz and a2 At In the same reference | |128| , deceleration is expected as the sys- 



^^ Lousto |l28| comments only this former part and not the acceleration boost taking place after 
the Zerilli potential peak. 
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tern is losing energy and momentum. This repulsive behaviour - before the Zerilli 
potential peak - is confirmed in the abstract of | |129| . 

Conversely, for Barack and Lousto 1 181 the radial component of the self- force is 
found to point inward (i.e., toward the black hole) throughout the entire plunge, in- 
dependently on the starting point ZuO- The work done by the self-force is considered 
positive, resulting in an increase of the energy parameter E throughout the plunge. 
To these results, it is attached a specific gauge choice (as opposed to the energy flux 
at infinity, which is gauge invariant) p8| . 

The upward or inward direction impressed on the particle by its mass and the 
emitted radiation and whether this direction is maintained throughout the plunge or 
part of it, is a fundamental if not the main one, feature to acquire in such analy- 
sis. Again here, the statements from the three papers might not be contradictory as 
Lousto |128[|129| describes motion in coordinate Schwarzschild time, and includes 
geodesic deviations. Conversely, Barack and Lousto 1 18] describe motion in proper 
time and apply only the self-force without geodesic deviations. Nevertheless, it is of 
interest to remark that the concept of repulsion resurfaces again solely in coordinate 
time as in the elder debate. The coefficient of the geodesic deviation coefficient ai 
changes sign during fall, while it does not occur to a2 |[7j| which remains negative 
throughout. 

For a particle starting from rest at a finite distance from the black hole, an analytic 
approximation of the self- force for the / > 2 modes (while for / = 0, 1 the solutions 
in |218| are mentioned) is given by Barack and Lousto |T8|: 



,rl=0 , ^rl=l , V J 15..2^' / Z.2 ■ 4M .\ 1 



1=2 



^;.,/=f;.r+f/;/r+i: -Y^«Vr +v-Mz2+^(^ ) ' ^^^^ 



where E is the orbital energy of the particle. The force has only negative and even 
powers of /, which makes the sum quickly convergent and provides an excellent 
approximation for the numerical evaluation of the first few lower multipoles. The 
derivation of such net expression for the self- force is not described in |18|p 



8.2 Regularisation parameters 

Regularisation parameters of the mode-sum method have been confirmed indepen- 
dently by different papers fTT,'22l and they are consistent with the results obtained 
by the application of the ^ function |18|. In radial fall, there is a regular gauge 
transformation between the de Bonder and Regge- Wheeler gauges pT| , and thus 
the regularisation parameters were also determined in the latter gauge |[5][T8J. The 
results are: 



^^ In the Rapid Communication [18] there are 7 citations of a yet unpublished material containing 
mathematical and numerical justifications of the results therein. The author acknowledges private 
communications by L. Barack. 
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(69) 


where / = 1 - 2M/z„ and z„ = 


-(£2_/)l/2. 





5. J Effect of radiation reaction on the waveforms during plunge 

Thanks to a suggestion of B. Whiting, prehminary indications were found (5). The 
waveforms shifts are of the order of (tens of) seconds for a particle sensibly radiating 
for few thousands of seconds, having started at rest from a finite distance between 
4M and 40M. The assumption used therein is energy-momentum balance (the en- 
ergy radiated to infinity and absorbed by the black hole is imposed to be equal to 
the energy change in the particle fall). This assumption is likely jeopardised by the 
lack of instantaneous energy conservation [1651. 

The correct alternative is the computation of the self- force and its continuous im- 
plementation all along the trajectory. It is thus mandatory to consider the application 
of the recently proposed self-consistent prescription [ |91| , unfortunately not yet part 
of the state of the art in terms of its application. 



9 Beyond the state of the art: the self-consistent prescription 



In | [9T| a rigorous derivation and application of the self-force equation is proposed. 
In the derivation, emphasis is put upon the de Bonder gauge and the consequences 
of relaxation, that is not enforcing this gauge. On one hand, the de Bonder gauge is 
imposed by the nature of the self-force, solely defined in this gauge. On the other 
hand, the relaxation of the de Bonder gauge stems from the need of departing from 
the background geodesic to find the self- force that causes such departure. Previous 
derivations were based on the assumption of deviations from geodesic motion ex- 
pected to be small; by consequence, the de Bonder gauge violation should likewise 
be small. Instead, the new derivation by Gralla and Wald is a rigorous, perturbative, 
result, obtained without the step of de Bonder gauge relaxation and containing the 
geodesic deviation terms. 

But Gralla and Wald go a step further f9r|. For the evolution of an orbit, 
rather than a first order perturbation equation containing geodesic deviation terms, 
they recommend a self-consistent approach. Such prescription basically affirms the 
greater accuracy of a first order perturbation expansion along a continuously cor- 
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reeled trajectory as opposed to a higher order perturbation expansion made on the 
background geodesicj^ Self-consistency bypasses the issue of relaxation, since at 
each integration step a new geodesic is founcf^ 

The 'classic' first order perturbative expansion for the motion of a small body 
determines that the first order metric perturbations satisfy: 

V^Vyhocp-2R\p^hys = -lenMuaUp^-^^, (70) 

where x^ = corresponds to a geodesic 7 of the background spacetime, and u^ is 
the tangent to 7. It is reminded that hg^p = hg^p — jhgg^p with h = hg^pg^^ . For a 
retarded solution to this equation (thus satisfying the de Bonder gauge condition) of 
the type: 

h"^p{x)=M£'^ (^G+p„,^, - \gapG^y^>p}j [xM^')] u"'uP'dt', (71) 

the first order in A deviation of the motion from 7 is expressed by (see |91 1 for the 
additional spin term): 

u^^jiuNpZ") = -Rp^s"u^zru^ - ig"P+u«uP)iVshf^-hphfsyu^. (72) 



Self-consistency prescribes that rather than using eqs. ( |7Q|71|72| ), it is instead 



preferable to apply the self- force coherently all along the trajectory: 

5(3) \^!i _^^(f)] ^r 

V^Vyh,,p-2R\p^hys = -l67lMua{t)up{t) ^ /^ ;77^ (73) 

uPVpu" = -{g«P+u"uP){Vshf^- \^ph';'sWu^, (74) 

/z-p(x) =m£J Tg+p^,^, - IsapG^^^p) [x,Zp(T')] u'^'u^'dx', (75) 



^^ The evolution of an orbit is lately getting the necessary concern. Pound and Poisson [161] apply 
osculating orbits to EMRI, but unfortunately their method is not applicable to plunge, for two 
reasons: the semi-latus rectum of the orbit, which decreases for radiation reaction, is smaller than 
a given quantity, considered as limit in their study case; the velocities and fields in the plunge are 
highly relativistic and their post-Newtonian expansion of the perturbing force becomes inaccurate. 
Non- applicability to plunge stands also for the work by Hinderer and Flanagan 1 103 1. 
^^ Indeed, it affirms that it is preferable to apply successively slV^ order expansion atxo and then at 
xi,X2, ...Xfn, rather then a 2^^ or higher order expansion at solely xq. It is evident, though, that self- 
consistency and perturbation order are decoupled concepts and that the former may be conceptually 
applicable to higher orders and more specifically, when, and if, a second order formalism will be 
available: in the same line of reasoning it would be preferable to apply successively a 2^^^ order 
expansion at xq and then ditxi,X2, ... Xm, rather then a 3^^ or higher order expansion at solely xq. 
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where this time u^{t) in eqs. ( 74|75 ), normaUsed in the background metric, refers 
to the s elf-consistent motion Zp(T), rather than to a background geodesic as in eqs. 
(|7l|72|); G^R(^/Qf is the retarded Green function, normaUsed with a factor of — 16;r 



|164J ; the symbol %^^ indicates the range of the integral being extended just short of 
the retarded time T^et, so that only the interior part of the light-cone is used. 



The geodesic deviations vanish in eq. (74), since self-consistency is imposed. 
Nevertheless, there might be situations where, for a whatever reason, the numerical 
implementation of the self-consistent prescription may be cumbersome. In this case, 
the addition of geodesic deviation terms might as in eq. ( [72| ) be necessary. 

If furthermore, motion is wished to be expressed in coordinate time, like in the 
pragmatic approach, it is useful to find a correspondence between the first order 
perturbation equation, containing geodesic deviations, eq. ([72|), and the terms of eq. 



p9\ . It is going to be shown that the a^ term of eqs. ([59] [62]) corresponds to the self- 
force term of eq. (|72|), when the latter is transferred to coordinate time. Referring to 



the time and radial components of the self-force, it is obtained from eqs. (46 47 ): 



Jt2 



^ ^self - 



- mFo u^ u^ ■- 



-mg 



//3 









ufu^ 



il^, (76) 






: F;,,f = -mr^^uPur-m/P (h^^.g - ^-h;s,p^ u-^u'-mW, {11) 



where 



Since: 



^''=(^^r;5-^^5r;/3l^''^^^'^^- 



(78) 



d _ dt d d^r _ dr d^t d^r f dt\ 

dx dx dt dT^ dt dT^ dt^ \dx ) ' 

after some computation, the term k^ disappears when the self-force is expressed in 
coordinate time: 



d^r 
m—-7r = m 
dt^ 



ny''''^s''[^h-,s-2^W,p 



v^v^ 



tu{t) 



mFn yP v^ — mg^^ 



^(5y;5~ 2^y5;(5 



vrv^ 



(79) 



Furthermore, a tedious computation shows that eq. ( [79] ) is nothing else than the a^ 
term of eq. ( 59 ) apart from the regularisations by mode-sum or Riemann-Hurwitz ^ 
function: 
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«6 ^ g'^ {h*p^,s - lh;s,p^ vV4w -S'^ {hl^;5 - \h;s,p) v'v'. (80) 

The recasting of the Riemann tensor term and of the left-hand side of eq. ( [72|) 
into coordinate time determines the relation to the a\Az and a2 At terms in eqs. |59] 



62 ). Finally, the geodesic deviation equatiorPleq. (72 ), is dealt with by Aoudia and 



Spallicci (tJ as difference between background and perturbed motions. 



10 Conclusions 

It has been shown that free fall is still the arena for a deeper comprehension of 
gravitation. Furthermore, it still generates acute observations like relativistic gliding 
for which the asymmetric oscillations of a quasirigid body slow down or accelerate 
its fall in a gravitational background | [95| . 

But is the problem of radial fall solved? Do we know the laws of motion of a star 
falling into a black hole, the relativistic modem version of the falling stone? A fair 
and objective answer leads to a moderate optimism. The general relativistic problem 
has had undeniable progress from 1997, but a careful analysis of the literature shows 
that some issues are either still partly open or simply not fully at hand, in terms of 
clear procedures, by means of which clearly cut answers are obtained. 

The remaining steps to be fulfilled for a satisfactory level of comprehension for 
radial fall of a small particle into a large mass (represented by a SD black hole) 
within the first perturbative order in m/M are divided in three groups. 

Investigations to be pursued before recurring to the self-consistent prescription: 

I Compute in proper time the first order perturbation equation with deviation 



terms, i.e. all terms of eq. (72). 



II Identify and compare expressions for renormalised, semi-renormalised and un- 
renormalised accelerations in the SD perturbed geometry; evaluate repulsive 
effects. The last century debate on repulsion has not addressed the effects of 
mass finitude (/ = 0, 1) and radiation reaction on the motion (/ > 2). The inclu- 
sion of perturbations blurs the issue further and renders it more complex. 

III Compare the waveform corrections obtained with energy-balance by Aoudia 
|[5| with those to be obtained by the integration in the equation of motion of the 
approximate analytic expression of the self- force of Barack and Lousto flSJ. 
The comparison is of interest only in case of concurring outcomes, as both 
procedures have a large degree of approximation. An other comparison | [189| 
could be made with the numerical results for a 10 : 1 mass ratio, as larger mass 
ratios imply the solution of the numerical two-scale problem, not yet at hand. 

IV Solve the initial conditions for particles with non-null or even relativistic initial 
velocities with and without large eccentricities. If the outcome of the previ- 
ous item is successful, coarsely identify possible signatures of radiation reac- 



See Levi-Civita fl2l][l22 , Ciufolini and Wheeler |37|, Ciufolini J36j. 
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tioij^ It is likely that a non-null initial velocity would act as an amplifier of 
the waveform shifts due to radiation reaction. In |132| , Lousto and Price show 
that the energy radiated by a non-null initial velocity may rise up to almost two 
orders of magnitude (see also fig. 3 in | [127| ). Incidentally, this line of work 
is confluent with the interests shown recently in the particle physics commu- 
nity |8, 183','190"201"2121 for head-on collisions and the associated radiation 
reaction |89,90|. Hopefully, progress in numerical relativity may lead to analy- 
sis of larger mass ratios and thereby test future results of perturbation theor}f^ 
in the range beyond 10~^. 

Investigations to be performed solely by the implementation of the self-consistent 
prescription: 

V Evaluate the trajectory by means of the Gralla and Wald self-consistent method 
||91 1. At each integration step, for a given number of modes: evaluation of the 
perturbation functions at the position of the particle; regularisation by mode- 
sum or ^ methods; correction of the geodesic and identification of the cell 
crossed by the particle; computation of the source term; reiteration of the above. 
Perturbations and self- force analysis in de Bonder's gauge are on-going |^19j 

VI Repeat steps 1-4 on the basis of the acquired self-consistent trajectory. An iter- 
ative scheme may be envisaged also in coordinate time. 

Investigations to be performed independently from the self-consistent prescrip- 
tion: 

VII Identify whether in a thought radial fall experiment, there is a physical observ- 
able independent of the gauge choice or else manifesting a recognisable effects 
in a given gauge, following Detweiler |56|. 

VIII Identify the domain of applicability of the self- force. Although self-consistency 
represents a closer description to perturbed motion, it is limited to cases where 
local deviations are small, since after all, it remains a perturbative approach. 
Thus, a quantitative identification of the domain of applicability of such pre- 
scription would be of interest to the community developing LISA templates. 
The adiabatic approximation can't be evoked to establish the limits of any self- 
force based analysis. Indeed, the hypothesis on the feeble magnitude of lo- 
cal deviations may be less constraining than the adiabatic hypothesis, for the 
absence of requirements on averaging. An explicit condition, referring to the 
geometry of the orbit and to the mass ratio and defining the domain of applica- 
bility of the self-force versus fully numerical approaches, is not yet available. 



^^ For Mino and Brink 1 148 1, the energy and momentum radiated are computed on the assumption 
that the small body falls in a dynamical time scale, with respect to proper time, well short of the 
radiation reaction time scale and therefore the gravitational radiation back-action on the orbit is 
considered negligible. The particle plunges on a geodesic trajectory, incidentally starting from a 
circular orbit, thus at zero initial radial velocity. 

^^ For head-on collisions. Price and Pullin have surprisingly shown the applicability of perturbation 
theory for the computation of the radiated energy and of the waveform for two equal black holes, 
starting at very small separation distances 1163), the so called close-limit approximation. 
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IX Investigate other gauges than Regge- Wheeler and other methods than the self- 
force for radial fall. A confirmation of the results by Barack, Lousto in (TSJ 
1 129) by an independent group is missing, though in |5| the necessary basic 
tools are developed, e.g. a numerical programme confirming the waveforms 



of |[T33l [T 3 9 



It would be also beneficial to carry out an analysis by the EOB 



method ||45 

The self- force community and even more the EMRI community are animated by 
different interests ranging from fundamental physics and theory, to numerical appli- 
cations, data analysis and astrophysics. To these variegated communities and gen- 
erally to physicists and astrophysicists, the capture of stars by supermassive black 
holes will bring a development comparable to the ancient markings made by the 
leaning Pisa tower, the Cambridge apple tree and the Einstein lift. 
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